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Abstract. This paper reviews current theoretical work on the various stages of accre- 
tion in protostars, and the relationship of these ideal stages to the spectral classes of 
0^ ' observed protostellar systems. I discuss scaling relationships that have been obtained 

Cn , for models of pre-stellar cores as they evolve by ambipolar diffusion toward a central 

singularity, and expectations for the dynamical evolution as the cores collapse radially 
to form a rotating disk. I summarize work that suggests accretion in T Tauri systems 
f"-., ' may be limited by ionization to disk surface layers, and implications for the variation 

CO . of Mjj with radius. Finally, I describe models for the asymptotic structure of free 

^\^ ' magnetocentrifugal winds, and show that the constant-density surfaces in these winds 

may be strongly coUimated even when the streamlines are not. 
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i^ ■ INTRODUCTION 

\Z( • The subject of this review is accretion in protostellar systems, which is properly 

c/3 . a hot topic in a figurative but not a literal sense. In the other accreting systems 

"^ ' discussed at this meeting, the observed high-energy emission often dominates the 

longer-wavelength fiux, and much of the current theoretical effort is devoted to 
^ ■ developing dynamical and radiative models of hot accretion flows that are capable 

of producing the observed spectral distributions. For low-mass protostars and T 
Tauri stars, on the other hand, the observed X-ray emission of 10^^ — 10'^^ erg s~^ 
[1-4] is much enhanced relative to that of main-sequence stars, but still constitutes 
just a small fraction (lO""^ — 10~^) of the sources' bolometric luminosities; it can be 
understood as thermal emission from T ~ 10^ — lO^K plasma in the stellar magne- 
tospheres and inner parts of the jets of these systems [5,6]. The bulk of the emission 
from young, low-mass stellar objects (YSOs) originates instead in a cooler accretion 
flow, where gas temperatures range from ~ lO'^ K (in the boundary layer/hot spot 
region), to a few xlO^K (in the inner, centrifugally-supported, accretion disk) to 
~ lOK (in the outer, radial parts of the flow). As will be described below, many 
of the current issues facing theorists interested in protostars have a quite different 
character from the issues under investigation in higher-energy systems. 

In this paper, I will begin by briefly summarizing the standard observational 



classification of low-mass YSOs, which assigns sources to one of four broad groups 
based on their spectral energy distributions. As Nuria Calvet describes in detail 
in her observational review in this volume, other properties of YSOs - such as the 
relative masses of their disks vs. envelopes, their accretion rates, their outflow 
momentum flux, etc. - also correlate with their spectral properties. The spectral 
classification is generally interpreted as an evolutionary sequence corresponding to 
successive stages of protostellar accretion [7-10]; this paper outlines current under- 
standing of the underlying theory. Since space constraints make a comprehensive 
review impossible (see instead the volume [11] and its successors), I will focus on 
some of the ways in which protostellar accretion is different from the other sorts 
of accretion discussed at this meeting. I will highlight some recent results in the 
theory of the early (i.e. largest-scale) stages of protostellar accretion, in which 
the flow is primarily radial. I will also describe recent work suggesting that even 
Shakura-Sunyaev-type "a" -models of the late, T-Tauri phase may be quite differ- 
ent from the standard (e.g. [12]) M^) = constant paradigm. In the last section, 
I will describe models of free, wide-angle YSO winds that are different from the 
confined, narrow jets found in high-energy systems (cf. the review of Blandford in 
this volume). 



OBSERVATIONAL CLASSES OF YOUNG STARS 



The spectral classification of YSOs introduced by Lada [13] and expanded by 
Andre, Ward-Thompson, & Barsony [14] divides them into classes O-III, as follows: 

Class 0: These are embedded objects whose spectral energy distributions (SEDs) 
peak in the sub-millimeter and show little or no IR, and no optical. They are 
believed to be the earliest protostars, since the masses of their envelopes 0.2 — 2M0 
exceed the masses of their disks + protostars [14-16]. From the relative numbers 
of these sources compared to other YSO, ages are probably a few xlO'^ years. 

Class I: These are embedded objects which are strong far-IR emitters, with a 
peak in the SED at 10 — 100/i. They are believed to correspond to late protostars 
or "pre- T-Tauri" stars since their envelope masses are typically 0.1 — 0.3Mq or less 
[15]. From the relative number of these sources, ages are estimated at a few xlO^ 
years (e.g. [17]). 

Class II: These are optically-revealed objects whose spectral energy distributions 
consist of a stellar photosphere plus significant "excess" IR and UV. These are 
interpreted as systems containing a star, an accreting disk, a boundary layer/hot 
spot, and a remnant envelope. Disk masses derived from mm continuum during 
this stage are 0.01-0.1 Mq [18] Ages of classical T Tauri stars (CTTS) (nearly 
coincident with this class) are 1 — 4 x 10^ yrs [19]. 

Class III: These sources present stellar blackbodies, with no sign of active ac- 
cretion; this class is essentially the same as the weak T Tauri stars (WTTS). 



STAGES OF PROTOSTELLAR ACCRETION 

Ambipolar Diffusion and Density Restructuring 

Molecular clouds in our Galaxy contain hundreds of thermal Jeans masses 
{Mj = c^(7r/G)^/^p^^/^); the fact that the observed star formation rate is orders of 
magnitude below the ratio of the gas mass to the gravitational collapse time is used 
to argue that molecular clouds are primarily supported by magnetic pressure and 
tension rather than thermal pressure (e.g. [7,20,21]). Under these circumstances, 
dynamical collapse to form a star occurs only when the ratio of mass to magnetic 
flux exceeds a certain critical value [22,23]; for the central flux tube the correspond- 
ing ratio of surface density to magnetic field is Tj/B = 1/{2tt\/G). A given column 
density for a cloud implies a certain critical magnetic field strength Bcru, and vice 
versa. Using a field strength of 20fiG, and a size scale of ~ O.lpc, the critical mass 
(for a spherical cloud) corresponds to ~ Mq, comparable to the minimum masses 
of dense cores where low mass stars are observed to be forming. 

As Mestel and Spitzer pointed out in the 1950s, the first stage in the protostellar 
accretion process under present Galactic conditions is expected to be the inward 
diffusion of neutrals relative to ions and the magnetic field, in order to create a 
supercritical core within a subcritical cloud. Thus the first accretion timescale 
is the ambipolar diffusion timescale tAo = L/vdrift, which can be estimated by 
equating the drag force on ions from inflowing neutrals with the magnetic force on 
ions; the result is [24] 

tAD^2x lOVrs [-J^zr) [-^) ■ (1) 

The ambipolar-diffusion/density restructuring stage has been studied by Nakano 
[25] and Lizano and Shu [26] in the quasistatic limit (applicable to early phases), 
and by Mouschovias and his collaborators [27-29] with a full dynamical treatment. 
For the latter numerical calculations, evolution is typically followed until the cen- 
tral density reaches 10^ cm~^, 1-2 orders of magnitude below the limit at which an 
opaque "protostar" would form. The same phase of evolution of initially subcritical 
cores has also recently been studied semi-analytically by modeling magnetic stresses 
in an angle-averaged sense [30]. In complementary studies, Tomisaka [31] has per- 
formed ideal-MHD (i.e. no ambipolar diffusion) simulations of the evolution of 
magnetized clouds which are initially supercritical, beginning from a cylindrically- 
stratified equilibrium. While the results of all these authors differ in detail depend- 
ing on their assumptions and parameter regime, the similarities in the results are 
more striking; I will concentrate on the latter. In both the initially-subcritical and 
the initially-supercritical simulations, the model clouds evolve towards flattened 
(disk-like), infalling structures in which the surface density S oc 1/R, the magnetic 
fleld strength B oc 1/R, and the midplane density pc oc l/R"^, where R is the cylin- 
drical radial coordinate. The peak inflow speed is 2-3 times the sound speed (which 
for T = lOK is c, = 0.2 km s'M. 



Schematically, these common properties can be understood from dimensional 
analysis as follows: First, marginally supercritical disks obey B = 2nG^''^Yl lo- 
cally, so that once a region of the cloud becomes supercritical it will maintain 
-B oc S because the ambipolar diffusion time is long compared to the dynam- 
ical time. Second, in local vertical equilibrium, an isothermal disk has central 
density pc = 7rGS^/(2c^) [32]; hence we can expect 5 oc S oc p^^ for the 
dynamically-evolving part of the flow. Finally, near-balance of the inward force 
of gravity (~ GM/R^ ~ GS) against the large-scale gradients of magnetic pres- 
sure and tension (~ B^/{pcR) -^ GT? / {p^R) — * c^/-R) and thermal pressure 
(~ pcCJ./{pcR) -^ cj./R) would imply S ~ cl/{GR), and hence B oc 1/R and 
Pc oc 1/ R^. Ideas along these lines have been developed with more rigor by Basu 
[33] (see also this volume), who constructs an approximate self-similar dynamical 
solution that matches well to his numerical solutions for this stage of accretion. 



Dynamical Collapse and Centrifugal Disk Formation 

What comes next? How does dynamical collapse proceed in the centrally- 
condensed structures, so as to appear as the observed Class and Class I objects? 
At this point, there is not the same convergence of theoretical models for the dy- 
namical collapse stage (in which speeds approach free-fall in the inner parts) as 
there is for previous restructuring stage (which involves speeds no greater than 
a few km s^^). Because the result of the previous stage yields a density profile 
p oc r~^, though, we may anticipate close resemblance to the family of (unmagne- 
tized) isothermal-sphere self-similar collapse models whose best-known representa- 
tives are the Larson [34]/Penston [35] and Shu [36] solutions. 

To gain insight into the dynamical collapse process, Galli & Shu [37] calculated 
solutions for collapsing clouds as perturbations with nonzero magnetic field strength 
about the inside-out singular isothermal sphere collapse. They found that, even 
with an initially uniform magnetic field and spherical density distribution, magnetic 
pinching forces lead to the formation of a growing equatorial "pseudodisk." A 
small-amplitude magnetic field, however, does not alter the accretion rate; this 
remains equal to the value 0.975c^/G of [36]. In an opposite limit, Li & Shu [38] 
demonstrate that an initially static disk with initial magnetic and surface density 
profiles -B oc S oc 1/R (i.e. the asymptotic scalings obtained from the evolution 
in the "restructuring" stage) undergoes an inside-out collapse. The accretion rate 
is given by M ^ (1 + Hq)cI/G, where 1 + ifo is the density overfactor allowed 
by magnetic support in the initial equilibrium compared to an unmagnetized state 
{Ho = 0). Although the calculation was performed in the Hq » 1 thin-disk limit, 
more realistic cases where Hq = 1 — 10 may be expected to show similar behavior for 
the dynamical collapse that follows the formation of a central density singularity. 

In the very innermost regions, the solution must approach free-fall, and the 
asymptotic collapse solution for a perfectly cold, unmagnetized disk with ini- 
tial surface density Ej oc R~^ is informative. For disk with initial profile 



Ej(i?) = Mj(i?)/(27ri?^) oc R ^ where Mi oc R, the resulting asymptotic scal- 
ings (for large time t and small radius R) are S oc t^^'^R^^'^, v oc t^'^i?^^'^, and 
M = const. = Mi{R)/tcoUapse{R) = 2^{nGy^^{^i{R)Rf/'^. Using the surface den- 
sity profile S ^ c1/{GR) of [29], one obtains M = Ucl/G, or about 2 x IQ-^Mg/yr 
for typical conditions. The associated collapse time ~ lO^yr is comparable to the 
observationally-estimated ages of embedded sources, giving some confidence that 
such dynamical collapse models may indeed be associated with the Class and 
Class I objects. 

As collapse proceeds, infalling material will finally create a centrifugal disk whose 
radius increases in time due to the larger specific angular momentum of collapsing 
material at later times. Previous analysis of the growth of the centrifugal disk 
(see [39-41]) may be modified in the future to account for the differences between 
magnetized and unmagnetized (exterior) radial collapse models. Although the the- 
ory of the dynamics in this accretion stage remains rather uncertain, observations 
seem to have identified at least one case - the source HL Tau - where a large-scale, 
radially-infalling "pseudodisk" imaged in gas [42] surrounds a much more compact 
(presumably centrifugal) disk imaged in dust [43,44]. 

Accretion in Protostellar Centrifugal Disks 

In the earlier stages of their evolution, protostellar disks are massive and self- 
gravity must be important to their evolution. A number of workers have inves- 
tigated analytically and numerically the growth, development, and saturation of 
large-scale, low-m gravitational instabilities (i.e. spiral density waves) in disks 
with protostellar properties (see e.g. [45] for recent calculations and references). 
Others have taken steps toward characterizing the action of gravitational torques 
by an effective "a" parameter (e.g. [46] and references therein); this approach is 
more fruitful for local instabilities than global disk modes. In addition, still others 
have discussed how the presence of binary companions (present in the majority of 
systems) or planets can induce accretion through externally-excited gravitational 
torques (e.g. [47,48]). At present, however, there is no generalized theory for how 
the processes of density-wave driven accretion and binary companion growth may 
develop in tandem or in competition out of self-gravitating disk perturbations. 
Significant future efforts will be needed before this stage of disk evolution is well- 
understood theoretically. 

In the latter stage of centrifugal disk accretion, on the other hand, disk self- 
gravity is less important, and we believe we understand the relevant processes 
better. At present, the leading contender for a physical mechanism to produce 
the "viscous" torques in accretion disks around compact objects is the saturated 
Balbus-Hawley instability (see the contributions of Balbus and Gammie in this 
volume). Does the same process work in the disks of YSOs as well? The issue 
here is whether, with the low temperature and high surface density conditions of 
protostellar disks, the gas is sufficiently ionized so that resistive dissipation does not 



suppress the instability. In the inner parts (i? ^ 0.1 AU) of protostellar disks, the 
central disk temperature exceeds lOOOK and collisional ionization of alkali metals 
provides sufficient conduction electrons. But in the low-temperature outer parts of 
these disks, some source of external ionization is needed. 

Recently, Gammie [49] argued that if cosmic rays alone provide the outer disk 
ionization, just the upper and lower surfaces to a depth of 100 g cm~^ will be 
sufficiently ionized to couple magnetic fields to the disk matter, and that the result 
would be a "layered model" where disk surfaces accrete and the equatorial portion 
is inert. In corresponding "alpha" models of disk accretion, the "active layer" 
surface density S^ is effectively constant with R, while the mass accretion rate M 
varies with R depending on the local opacity. The result is a relatively low value 
of accretion 






M = 1.8 X lO^^Meyr"^ (2) 

°-^ VO.Ol/ VlOOgcm-3; ^ ' 

in the innermost region, and the possibility of mass accumulation in the outer disk if 
the infall rate exceeds the carrying capacity of the surface layers. The low predicted 
value of M in the "layered" model appears consistent with recent spectral modeling 
of observed classical T Tauri disks (see the chapter by Calvet). If the difference 
between the infall rate and accretion rate is large enough, the mass buildup in 
the outer disk could periodically be purged via gravitational instabilities; such a 
mechanism could potentially be responsible for the observed FU Ori outbursts (see 
the chapter of Calvet for discussion of these phenomena) . 

Following up on these ideas, Glassgold, Najita, & Igea [50] have shown that since 
YSOs are strong X-ray emitters, the surface layers of their accretion disks will be 
ionized to a thickness comparable to Gammie's estimate, even if the cosmic ray flux 
is inhibited by strong, magnetized T Tauri winds. Further work on these sorts of 
"layered" models will be needed to determine whether both the observed low mean 
accretion rate of T Tauri stars and their high-accretion-rate outbursts events can 
be explained within a unified framework. 

PROTOSTELLAR WINDS 

My final topic is the nature winds from protostellar systems, and their con- 
nection to the narrow Herbig-Haro jets and broader massive molecular outfiows 
which, observationally, seem to be an inescapable side-effect of protostellar ac- 
cretion. Optically-visible Herbig-Haro objects are believed to represent shocked 
regions in dense, high-velocity gas that emerges from the very innermost part of 
protostellar accretion disks (see contributions by Frank and Bally in this volume). 
Although there is general agreement that the surrounding molecular outfiows con- 
sist of ISM material that has been mobilized by an interaction with a wind, there 
is equally general disagreement about the nature of that primary wind, and its 
relationship to observed H-H jets. 



In most current models, the primary wind is believed to be magnetocentrifugally 
driven, from the inner disk edge where the accretion flow interacts with the stellar 
magnetosphere (see [51] and the review of Stone in this volume), and/or from the 
surface of the disk at larger radii (see the review [52]). The controversial issue is 
whether the primary wind consists solely of an isolated, highly coUimated flow - 
seen as the ionized, optical jet - or whether there is also a lower-density, wider- 
angle, neutral wind surrounding the observed jet. 

In the first "isolated jet" picture, molecular outflows are essentially the wings of 
bow shocks in the jet head or beam (see e.g. [53]); hence the gas motions transverse 
to the jet are driven by pressure gradients. One difficulty with such models is that 
the cooling under interstellar conditions may be too strong to allow a very broad 
bow shock to form (e.g. [54,55]). Another potential difficulty is that, depending 
on the variation of magnetic field strength with radius at the base of the flow, 
it may not even be possible to drive a wind which is both strongly collimated 
and fast [56]. A related difficulty is that since magnetocentrifugal wind models 
typically have internal Alfven speeds which approach (or even exceed) their flow 
speeds, they may not be able to remain collimated as they propagate into the 
low-pressure ambient medium - where the sound speed is 1-3 orders of magnitude 
smaller than the wind speed. Adam Frank's contribution to this volume describes 
ongoing numerical studies of jet/cloud interactions which address these and other 
issues. 

In the second, "wide-angle wind" picture, molecular outflows are interpreted as 
the pileup from a momentum-driven "snowplow" [57]. A wide-angle wind would 
occur if a magnetocentrifugal outflow is unable to self-collimate, or becomes decolli- 
mated as it emerges from a dense core into a medium of low ambient pressure. Such 
a wind would expand laterally to flll the whole solid angle between the equator and 
whatever natural boundary exists near the pole. Because cold magnetocentrifugal 
outflows only occur on fleld lines having angles less than 60° with respect to the 
equator, it is natural to expect that there would remain some unloaded axial flelds 
fllling the space interior to the wind towards the pole, which originate in the disk at 
angles greater than 60° with respect to the equator. In particular, when the wind 
is formed by tearing open a magnetosphere [51,59], the axial magnetic flux interior 
to the wind will be comparable to the poloidal magnetic flux in the wind itself. If 
the wind fleld lines are instead open fleld lines carried in with the accretion flow, 
then the ratio of interior magnetic flux to wind magnetic flux would depend on the 
distribution of the fleld in the disk. 

The asymptotic structure of a wide-angle MHD wind, for the case of a flow along 
opened magnetospheric fleld lines emerging nearly isotropically from the corotation 
region (the "x-wind"), has been analyzed by Shu, Najita, Ostriker, & Shang [60]. 
They show that, while the streamlines collimate only slowly with distance, the 
density structure in the wind becomes nearly cylindrically stratifled. Briefly, the 
reasons for these results can be explained as follows: In the asymptotic regime, the 
largest term in the internal stress is that associated with the toroidal magnetic fleld 
B^p. In order to have a locally force-free conflguration, d{B^R)/d6 ~ must hold. 



where R = r sin 6 is the cyhndrical distance from the axis, and 6 is the polar angle. 
This implies 

B^ = -C{r)/R (3) 

for some (positive) function C{r). From angular momentum conservation and field 
freezing, we have 

-^ = and Bp = PpVp, (4) 

Bp Vp 

where Vp and Bp are respectively the poloidal flow speed and poloidal magnetic 
field, Q is the angular rotation rate of the footpoint of the local field line, p is the 
density, and f3 represents the ratio of magnetic fiux to mass fiux (which remains 
constant on any streamline). From these equations, we find 



Cjr) 
Qf3R^ 






Using V ■ Bp = 0, it can be shown that C is a slowly- varying function of r. Thus, if 
(3 and Q vary little between streamlines (as is true in the x-wind solutions of Najita 
and Shu [61]), then the density in the wind will vary approximately as p oc R~^ at 
large distances. Note that while the scaling p oc r~^ arises kinematically from the 
continuity equation for a constant-speed radial fiow, the scaling p oc (sin 9)"^ comes 
about because hoop stresses redistribute the streamlines in 6 until the magnetic field 
is force-free (cf. eq. 3). The actual value of the function C{r) may be evaluated by 
balancing the pressure at the inner edge of the wind -B^/(87r) against the pressure 
provided by interior axial fields, -Bp j„(/(87r), and applying the constraint of mass 
outfiow conservation. The value of the poloidal speed on any streamline can be 
obtained from energy conservation; in terms of the streamline's footpoint radius 
Ro and Alfven radius Ra, Vp = n{2R\ - SR^y/^. 

For any variation of Q, (3, Rq, and Ra with magnetic fiux $ (this variation must be 
derived self-consistently for the wind inside the fast-MHD surface), the foregoing 
considerations can be applied to find the asymptotic structure of a "free" MHD 
wind - whether the wind originates in a narrow or broad region of the disk. It can 
be shown, for example, that disk winds which leave the disk with p oc Rq*^ and B oc 
Rq also have slowly-coUimating (i.e. pointing close to radially) streamlines 

but nearly cylindrically-stratified density structure near the pole, at large distances 
from the disk. Figure 1 shows an example of the asymptotic structure of a wind 
in which all streamlines originate at the same radius Rq and have Ra = 2i?0j and 
for which the interior axial magnetic fiux along the pole equals the wind poloidal 
magnetic fiux. 

Wide-angle "snowplow" models with winds of this sort have been successful at 
explaining many general characteristics of outfiows [62] , as well as particular outfiow 
systems [63]. 



Further work on both observational and theoretical fronts is still needed to iden- 
tify which set of conditions leads to a wide-angle wind, and which to an isolated jet, 
or perhaps to conclude instead that alternative models are warranted to represent 
the primary winds from protostars. 

I am grateful to Charles Gammie and Jim Stone for their helpful comments on 
this manuscript. 



REFERENCES 

1. Casanova, S., Montmerle, T., Feigelson, E. D., & Andre, P. 1995, ApJ, 439, 752 

2. Koyama, K., Hamaguchi, K., Ueno, S., Kobayashi, N., & Feigelson, E. D. 1996, 
PASJ, 48, L87 

3. Montmerle, R. Feigelson, E., Bouvier, J., & Andre, P. 1993, in Protostars and Planets 
III, E. Levy & J. Lunine eds, U. Arizona Press, p. 689 

4. Neuhauser, R., Sterzik, M. F., Schmitt, J. H. M. M., Wichmann, R., & Krautter, J. 
1995, A&A, 297, 391 

5. Feigelson, E. D., Gianipappa, M. S., & Vrba, F. J. 1991, in The Sun in Time, C. 
Sonnet & M. Giampappa, eds., (Tucson:Univ. Arizona Press), p. 658 

6. Feigelson, E. D. 1997, in Star Clusters and Age Indicators, S. Sciortino &; R. Pal- 
llavicini, eds. 

7. Shu, F. H., Adams, F. C, Lizano, S. 1987, ARAA, 25, 23 

8. Adams, F. C, Lada, C. J., & Shu, F. H. 1987, ApJ, 312, 788 



ISl 



1000 
800 
600 
400 
200 




- 


' / / / '- 


- 


\\\ / \ 


1 / / '/ / + 


['//// '^ 


: u / /'i / /^ 


/// / / /-^^^ '- 


\\jl//\ ^^-^ \ 


V y / ^'^'^ * 


|//y^^^l t 


tr^ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



10* 
8000 
6000 
4000 
2000 






200 400 600 800 1000 

R/R„ 



2000 4000 6000 8000 ^^^ 

R/R„ 



FIGURE 1. Asymptotic magnetocentrifugal wind, on two scales. Strcamlines/ficld lines are 
shown with solid curves at equal increments of mass loss/magnetic flux. Density contours are 
shown with dashed curves, for log{p/ po) = —4, —5, —6, —7 (left), and log{p/ po) = —6, —7, —8, —9 
{right); po = M^/{iTTRln). 



9. Lada, C. J. 1995, in Molecular Clouds and Star Formation, C. Yuan &: J. You eds, 
(Singapore: World Scientific), p. 1 

10. Shu, F. H. 1995, in Molecular Clouds and Star Formation, C. Yuan & J. You eds, 
(Singapore: World Scientific), p. 97 

11. Levy, E. & Lunine, J. (eds) 1993, Protostars and Planets III, (Tucson:U. Arizona 
Press) 

12. Pringle, J. 1981, Ann. Rev. Ast. Ap., 19, 137 

13. Lada, C. J. 1987, in lAU Symposium No. 115: Star Forming Regions, M. Peimbert 
& J. Jugaku, eds. (Dordrecht:Reidel), p. 1 

14. Andre, P., Ward-Thompson, D., k Barsony, M. 1993, ApJ, 406, 122 

15. Andre, P., & Montmerle, T. 1994, ApJ, 420, 837 

16. Bontemps, S., Andre, P., Tereby, S., & Cabrit, S. 1996, A&A, 311, 858 

17. Wilking, B. A., Lada, C. J., & Young, E. T. 1989, ApJ, 340, 823 

18. Beckwith, S. V. W., Sargent, A. I., Chini, R., & Gusten, R. 1990, AJ, 99, 924 

19. Cohen, M., & Kuhi, L. V. 1979, ApJ Supp, 41, 743 

20. McKee, C. F., Zweibel, E. G., Goodman, A. A., & Heiles, C. 1993, in Protostars and 
Planets III, ed. E. Levy & J. Lunine (Tucson: University of Arizona Press), p. 327 

21. Ostriker, E. C. 1997, in Star Formation, Near and Far, eds. S. S. Holt and L. G. 
Mundy, (Woodbury, NY:AIP Press), p. 51 

22. Mouschovias, T., & Spitzer, L. 1976, ApJ, 210, 326 

23. Tomisaka, K., Ikeuchi, S., &: Nakamura, T. 1988, ApJ, 335, 239 

24. Mouschovias, T. 1989, in The Physics and Chemistry of Interstellar Molecular 
Clouds, G. Winnewisser & J. T. Armstrong, eds (Berlin:Springer), p. 297 

25. Nakano, T. 1979, PASJ, 31,697 

26. Lizano, S., & Shu, F. 1989, ApJ. 342, 834 

27. Fiedler, R. A., & Mouschovias, T. 1993, ApJ, 415, 680 

28. Ciolek, G. E., & Mouschovias, T. 1994, ApJ, 425, 142 

29. Basu, S., & Mouschovias, T. 1994, ApJ, 432, 720 

30. Safier, P., McKee, C, & Stabler, S. 1997, ApJ, 485, 660 

31. Tomisaka, K. 1996, PASJ, 48, 701 

32. Spitzer, L. 1978, Physical Processes in the Interstellar Medium (New York:Wiley) 

33. Basu, S. 1997, ApJ, 485, 240 

34. Larson, R. B. 1969, MNRAS, 145, 271 

35. Pension, M. V. 1969, MNRAS, 144, 425 

36. Shu, F. H. 1977, ApJ, 214, 488 

37. Gain, D., & Shu, F. H. 1993, ApJ, 417, 220 

38. Li, Z.-Y., & Shu, F. H. 1997, ApJ, 475, 237 

39. Cassen, P., & Moosman, A. 1981, Icarus, 48, 353 

40. Cassen, P., & Summers, A. 1983, Icarus, 53, 26 

41. Tereby, S., Shu, F. H., & Cassen, P. C. 1984, ApJ, 286, 529 

42. Hayashi, M., Ohashi, N., & Miyama, S. M. 1993, ApJ, 418, L71 

43. Mundy, L. G. et al. 1996, ApJ, 464, L169 

44. Looney, L. W., Mundy, L. G., k Welch, W. J. 1997, preprint 

45. Laughlin, G., Korchagin, V., & Adams, F. C. 1997, ApJ, 477, 410 

46. Gammie, C. 1997, in preparation 



47. Lin, D. N. C, & Papaloizou, J. C. B. 1993, in Protostars and Planets III, E. Levy 
& J. Lunine eds, U. Arizona Press, p. 749 

48. Ostriker, E. C, Shu, F. H., & Adams, F. C. 1992, ApJ, 399, 192 

49. Gammie, C. F. 1996, ApJ, 457, 355 

50. Glassgold, A. E., Najita, J., & Igea, J. 1997, ApJ, 480, 344 

51. Shu, F., Najita, J., Ostriker, E., Wilkin, F., Ruden, S., Lizano, S., 1994, ApJ, 429, 
781 

52. Konigl, A., & Ruden, S.P. 1993, in Protostars and Planets III, , E. Levy & J. Lunine 
eds, U. Arizona Press, p. 641 

53. Chernin, L., Masson, C, Pino, E.M., & Benz, W. 1994, ApJ, 426, 204 

54. Blondin, J.M., Fryxeh, B.A., & K6nigl, A. 1990, ApJ, 360, 370 

55. Stone, J.M, & Norman, M.L. 1993, ApJ, 413, 198 

56. Ostriker, E. 1997, ApJ, 486, 291 

57. Shu, F., Ruden, S.P., Lada, C, & Lizano, S. 1991, ApJ, 370, L31 

58. Blandford, R. D. & Payne, D. G. 1982, MNRAS, 199, 883 

59. Ostriker, E. C. & Shu, F. H. 1995, ApJ, 447,813 

60. Shu, F., Najita, J., Ostriker, E., & Shang, H. 1995, ApJ, 495, L155 

61. Najita, J., & Shu, F. 1994, ApJ, 429, 808 

62. Li, Z.Y., & Shu, F.H. 1996, ApJ, 472, 211 

63. Nagar, N. M., Vogel, S. N., Stone, J. M., & Ostriker, E. C. 1997, ApJ (Letters), 482, 
L195 



